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ABSTRACT 

The total mass of clusters of galaxies is a key parameter to study massive halos. It relates to numerous gravitational and baryonic 
processes at play in the framework of large scale structure formation, thus rendering its determination important but challenging. 
From a sample of the 1 1 X-ray bright clusters selected from the EXCPRES sample, we investigate the optical and X-ray properties 
of clusters with respect to their total mass derived from weak gravitational Iensing. From multi-color wide field imaging obtained 
with MegaCam at CFHT, we derive the shear profile of each individual cluster of galaxies. We perform a careful investigation of all 
systematic sources related to the weak Iensing mass determination. The weak Iensing masses are then compared to the X-ray masses 
obtained from the analysis of XMM-A'ewton observations and assuming hydrostatic equilibrium. We find a good agreement between 
the two mass proxies although a few outliers with either perturbed morphology or poor quality data prevent to derive robust mass 
estimates. The weak Iensing mass is also correlated with the optical richness and the total optical luminosity, as well as with the X-ray 
luminosity, to provide scaling relations within the redshift range 0.4 < z < 0.6. These relations are in good agreement with previous 
works at lower redshifts. For the Lx-M relation we combine our sample with two other cluster and group samples from the literature, 
thus covering two decades in mass and X-ray luminosity, with a regular and coherent correlation between the two physical quantities. 

Key words. Gravitational Iensing: weak - X-rays: galaxies: clusters - Cosmology: observations - Cosmology: dark matter - Galaxies 
: clusters : general - Galaxies : clusters : individual 



1. Introduction 

Clusters of galaxies have been attracting considerable interest 
for their cosmological applications since a few decades. The ob- 
served cluster abundance, converted into a mass function, can be 
confronted to cosmological numerical simulations to put tight 
constraints on cosmological parameters like the amplitude of 
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practical determination of the masses themselves. In the sim- 
plest model of structure formation and evolution which involves 
a pure gravitational collapse of dark matter halos, groups and 
clusters are expected to form a population of self-similar objects 
characterized by simple relations linking the total mass to other 
physical quantities ( |Kaiser||1986^ . These scaling relations rep- 
resent an efficient way to convert simple observables to masses 
for large samples of objects. However, the link between the to- 



* Based on observations obtained with MegaPrime/MegaCam, a joint 
project of CFHT and CEA/DAPNIA, at the Canada-France-Hawaii 
Telescope (CFHT) which is operated by the National Research Council 
(NRC) of Canada, the Institut National des Science de I'Univers of the 
Centre National de la Recherche Scientifique (CNRS) of France, and 
the University of Hawaii. 



tal mass of an object and its baryonic tracers is not fully under- 
stood. Several processes break the self-similarity and introduce a 
scatter around the theoretical relation that needs to be accounted 
for when determining the mass function (see |Voit||2005| for a 
review). Hydrodynamic simulations offers a way to complete 
this model. They can define more realistic relations between 



the cluster total mass and the X-ray observables (e.g. Kravtsov 
|et al . 2006; Nagai et al. 2007 ). But the physics of the baryons 
within these simulations is not well constrained. Therefore, it is 
of prime importance to characterize and calibrate scaling rela- 
tions with observational results based on accurate direct mass 
measurements. Such empirical relations can also improve our 
understanding of several non-gravitational processes disturbing 
the evolution of the large-scale structures. 

Among the several ways to derive the total mass of a cluster 
of galaxies, an efficient one is to use the X-ray emission of the 
intra-cluster gas. The distribution in density and temperature can 
be measured and, assuming the hydrostatic equilibrium of the 
cluster, it is possible to derive the dark matter potential shape. 
This kind of approach has been applied successfully on many 
X-ray observations ( Pointecouteau et al. ' 2005[ |Vikhlinin et al.| 
2006 ; Buote et al.|2007| l. However, masses derived from the X- 
ray measurements suffer systematics from non-relaxed objects 
or from non-thermal pressure support, and tends to underesti- 
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mate the true total mass of clusters by~ 10 - 15% as seen in 
numerical simulations (Kay "2004"; 'Nagai et al.|2007[ Lau et al. 



|2009 ; Meneghetti et al. 20 10) and suggested by observational 
results (e.g. |Mahdavi et al.|20 08 ). 

Gravitational lensing is another way to derive the total pro- 
jected mass, without any assumption about the dynamical state 
of the cluster The strong lensing regime that occurs in the cen- 
ter of some clusters leads to elongated arcs and multiple images 
of background galaxies which probe the central mass distribu- 
tion. In parallel, the weak lensing regime leads to estimates of 
the total mass through the statistical distortion of background 
galaxies up to large projected distances from the cluster center 
(e.g. Bartelmann & Schneider,2001, for a review). However, this 
method gives reliable results only for massive clusters in which 
the weak lensing signal is detectable with high enough accuracy. 
It is also limited to intermediate redshift clusters for which the 
distant background galaxies are far enough behind the cluster 
Gravitational lensing masses are also subject to systematics such 
as the determination of the redshift distribution of the sources or 
projection effects of matter structures near the cluster or along 



the hne-of- sight ( |Metzleretal.|2001[[Hoekstra|2003]|Meneghetti 
let al.|2010| l. 

Because all these mass estimators have their own weaknesses 
and limits and because they fail to properly recover the total 
mass in specific cases, there has been attempts to provide joint 
analysis of the mass profiles, combining strong and weak lens- 
ing (Kneib et al. 2003; Bra dac et al.| 2005; Merten et al. 2009; 
Oguri et al.|201 1,) or lensing and X-ray measurements ( Mahdavi, 
et al.||2007[ [Morandi et"aL][20T0) l. But the difficulty to get data 



of similar "quality" to combine them efficiently has up to now 
limited the use of such analysis. 

In order to better quantify the scaling relations of clusters 
with weak lensing masses, we present in this paper the study 
of a sample of 1 1 intermediate redshift clusters which are part 
of the EXCPRES sample. These clusters have been observed in 
X-ray with XMM-Newton (Arnaud et al. in preparation) and in 
wide field optical imaging with Megacam at the Canada-France- 
Hawaii Telescope ( Foex et al.|2012 hereafter Paper 2). The pa- 
per is organized as follows. In section 2 we present the sample of 
1 1 clusters selected for this study. Section 3 is dedicated to the 
lensing methodology with special attention to the shape mea- 
surements of the background galaxies and the normalization of 
the lensing signal. Lensing masses are determined from the shear 
signal in section 4 where we also make the comparison with the 
X-ray masses. Scaling relations are built and discussed in section 
5 and some conclusions are drawn in section 6. In all the pa- 
per, we assume a standard A-CDM cosmology with Q.m = 0.3, 
Qa = 0.7 and Hq = 70 hjo kms^' Mpc^'. 

2. The cluster sample description 

The 11 clusters used in the present study are a sub-sample 
of the EXCPRES sample {Evolution of X-ray galaxy Cluster 
Properties in a REpresentative Sample, Arnaud et al. in 
preparation). This is a representative X-ray selected sample of 
clusters in the redshift range 0.4 < z < 0.6, compiled from 
published X-ray selected samples at the end of 2004 (e.g., 
EMSS, NORAS, REFLEX, B-SHARC, S-SHARC, 160SD and 
WARPS-I samples). EXCPRES is designed to study the evo- 
lution of the X-ray properties of clusters. Objects are selected 
only on the basis of their X-ray luminosities, as measured in 
the aforementioned surveys, and their redshift. The EXCPRES 
clusters are selected in the Lx - z plane around a median 
redshift of z = 0.5, so that there are about 5 clusters in roughly 



equal logarithmically-spaced luminosity bins (factor of ~ 2 
variation between bins). No selection is added with respect to 
the dynamical state, ensuring that no bias is introduced by the 
sample construction process, making it ideal for the study of 
scaling relations and their evolution. The final sample contains 
20 clusters. 

The main issue to deal with such a sample of clusters is the 
selection bias. The 3 clusters from the surveys with a strict 
X-ray surface brightness limit (B-SHARC and REFLEX) are 
at 1-2 times the survey limit, i.e. in the regime where the bias 
is not negligible. For the other clusters from the EMSS and 
NORAS surveys, that have no defined flux limits, it is difficult 
to say wether there is a bias. Note however that NORAS 
clusters are below the 3e"'^ flux Umit where the survey is 
50% complete, so likely introducing a bias. Therefore, there is 
most likely a Malmquist bias, that in principle should be taken 
into account for the study of the Lx - M relation. However it 
is basically impossible to estimate precisely this bias due to 
the variety of cluster surveys from which we constructed the 
EXCPRES sample. Note also that both the Lx - M relations 
from]Rykoff et al.|(|2008|) and|Leauffiaud et al.| ( [^TO| ) (compared 
and combined with our results in section 5.3) are not corrected 
for the Malmquist bias. This issue represents a limitation for 
precise studies of the Lx - M scaling relation. 
However for the Mx - M^l mass relation, the original main 
goal of this paper, the Malmquist bias effect should not be 
important. It could play a role via possible residual bias in 
dynamical state: X-ray surveys preferentially select highly 
peaked cool core (thus relaxed) clusters close to the survey flux 
limit because the are over-luminous for their mass. Thus the 
sample may contain more relaxed clusters than the underlying 
population. This may affect the Mx - Mwl relation as one expect 
better agreement for relaxed objects. However in view of the er- 
ror bars and dispersion (see Figure|6]l, this is not the major woiTy. 

Note that EXCPRES is deliberately built in a similar way 
to the local representative sample REXCESS (Bohrmger et al. 
|2007 1; so the scaling and structural properties of REXCESS 
dCroston et al . 2008; Pratt et al. 2009; Arn aud et al 2010) will 
be used as a local reference for EXCPRES. 

Only clusters with an X-ray luminosity > 5 x lO"*"* erg/s 
in the [0.5-2.0] keV band within the detection radius were se- 
lected for an optical follow-up. This threshold in Lx ensures that 
the total mass of each cluster is high enough to provide signifi- 
cant weak lensing signal. These measured XMM-A^evvfo« lumi- 
nosities were converted into the [0.1-2.4] keV energy band for 
furffier analysis and comparison with other works (see Sec. 5.3). 
They are reported Table [T] together with global temperatures de- 
rived within apertures maximizing the signal-to-noise ratio for 
the XMM-Newton spectroscopic analysis. In the framework of 
the EXCPRES project, the weak lensing measurement of these 
11 clusters allows a one-to-one comparison of the X-ray and 
weak lensing mass proxies to assess the reliability of the total 
mass estimate from X-rays data only. 

3. Weak lensing 

Gravitational lensing distorts the intrinsic image of the back- 
ground galaxies. This is usually quantified by the change in the 
shape parameters of the sources. But the intrinsic shape is con- 
volved by the effect of the atmosphere and the instrumental dis- 
tortion before it is measured on CCD images. The PSF correc- 
tion is therefore a critical step in the weak lensing analysis and 
it has to be done with great care. This is usually validated thanks 
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Table 1. General properties of the clusters. Columns: (1) Cluster Name. (2,3) Equatorial coordinates of the X-ray peak. (4) Redshift. 
(5) X-ray luminosity in the XMM-Newton [0.1-2.4] keV energy band within the detection radius, Rjef. (5) Detection radius. (6) 
Global spectroscopic XMM-Newton temperature, non core excised. (7) Detection of strong lensing in the CFHT data. 



Cluster 


RA 
(J2000) 


Dec 
(J2000) 


z 


Lx[0.1,2.4]keV 

(h-^ lO^'-'ergs-') 


^det 

(arcmin) 


kT strong lens 
(keV) (Y/N) 


MS 0015.9+1609 


00'' 18" 


33.26' 


+ 16°26'12.9" 


0.541 


16.3 H 


t 0.1 


5.1 


9.0: 


tO.2 


N 


MS 0451.6-0305 


04'' 54" 


10.85' 


-03°00'57.0" 


0.537 


15.2 H 


tO.l 


3.7 


8.8: 


tO.3 


Y 


RXC J0856.1+3756 


08''56" 


12.69' 


+37°56'15.0" 


0.411 


4.25 H 


tO.l 


3.8 


6.9: 


t04 


N 


RX J0943.0+4659 


09*42" 


56.60' 


+46°59'22.0" 


0.407 


4.56 H 


tO.l 


4.1 


54: 


tO.2 


N 


RXC J1003.0+3254 


10''03" 


04.62' 


+32°53'40.6" 


0.416 


2.80 H 


t 0.2 


3.5 


3.6: 


tO.2 


N 


RXJl 120.1+4318 


11''20" 


07.47' 


+43°18'06.0" 


0.612 


5.82 H 


tO.3 


2.1 


5.2: 


tO.3 


N 


RXC J1206.2-0848 


12''06" 


12.13' 


-08°48'03.6" 


0.441 


20.7 H 


tO.l 


64 


9.5: 


tO.2 


Y 


MS 1241.5+1710 


12'' 44" 


01.46' 


+ 16°53'43.9" 


0.549 


6.86 H 


tO.2 


3.5 


4.6: 


tO.l 


N 


RX J1347.5-1144 


13'' 47" 


32.00' 


-11°45'42.0" 


0.451 


41.8 H 


tO.l 


4.9 


11.4 


±0.2 


Y 


MS 1621.5+2640 


16'' 23" 


35.16' 


+26°34'28.2" 


0.426 


5.09 H 


tO.9 


4.2 


5.8: 


tO.7 


Y 


RX J2228.5+2036 


22'' 28" 


33.73' 


+20°37'15.9" 


0.412 


13.3 H 


t 0.1 


8.8 


7.8: 


tO.3 


N 



to realistic simulations of data. We describe below the details of 
how we implemented our weak lensing pipeline up to the mea- 
sure of the total mass of the clusters. 



3. 1 . Objects selection 

The very first step of any weak lensing analysis consists in de- 
tecting and sorting galaxies and stars. This is explained in details 
in Paper 1, a methodology that is summarized as follows. Using 
SExTRACTOR in dual mode, we detect the objects on the x^,., im- 
age while astrometric and photometric parameters are measured 
on each individual g', r', /' and z' images. Stars, galaxies and 
false detections are sorted according to several criteria: posi- 
tion in the magnitude/central flux diagram which displays the 
star branch, size with respect to the size of the PSF, and stellar- 
ity index, according to the CLASS_STAR parameter. Objects in 
masked area are also removed. After this selection step, galaxy 
densities are of the order of 30 arcmin"^, and the completeness 
magnitude ranges from 24.5 to 25 in the r' band, depending on 
the data quality in the different cluster fields. 

In order to correct for the PSF smearing, we use the Im2shape 
software (Bridle et al. 2002 1, which has already demonstrated its 
potential to recover the intrinsic shape of the galaxies fCypriano' 
[eFTl. 2004; Bardeau et al. 2005, 2007 ; ILi mousin et al. 2007a b). 
Starting with a given model for the shape of each object, the 



code convolves it with the local estimate of the PSF. For sim- 
plicity, both the PSF and the object are modeled with a single 
elliptical gaussian profile. Exploring the space parameters with 
a MCMC sampler, the most likely model is found by minimiz- 
ing the residuals, the whole distribution being used to estimate 
robust statistical eiTors on each parameter. 

We first determine the local PSF by looking at the shape 
of the stars. The resulting PSF field is cleaned and smoothed 
by looking at the 10 nearest stars at each point and removing 
those that differ by more than 1 .5cr from the local average shape. 
This ensures to avoid any jump in the PSF pattern, e.g. near 
masked area. We also exclude objects having an ellipticity larger 
than 0.2, thus reducing the contamination of the stars catalog by 
false detections and faint galaxies. The PSF map over the whole 
field of view is then obtained by averaging the ellipticities of 
the 5 nearest stars at each galaxy position. We checked that our 
Im2shape implementation can recover point-like objects by ap- 
plying this PSF correction to each star Figure |2] shows the re- 
sults of the shape parameters measurements for these stars: the 



size distribution is dominated by point-sources and the orienta- 
tion is more uniformly distributed after the PSF coiTection. It 
is also shown Figure ?? where the regular pattern of the stars 
ellipticities due to the PSF (top row) disappear after the decon- 
volution. 



3.2. Validation of the shape measurements with simulated 
data 

To check the performance of our Im2shape implementation, we 
tested our pipeline on the simulated fields provided by the STEP 
collaboration (the Shear Testing Programme, Heymans et al. 
2006). STEPl simulations were built in order to compare the 



accuracy of different lensing pipelines for ground-based images. 
They provide a set of images reproducing CFHT-like observa- 
tions with representative densities of stars and galaxies. These 
simulations combine several PSF models and shear intensities. 
We ran our lensing pipeline on each corresponding data set, com- 
pare our results to the input configurations and derive an average 
shear calibration bias and PSF residual. 

Our implementation of Im2shape differs slightly from the one 
proposed by S. Bridle in the STEPl simulations: for simplic- 
ity and efficiency we fitted both the PSF and the galaxy shapes 
by a single elliptical gaussian instead of 2 concentric ones. This 
avoids to deal with 2 different values of the galaxies ellipticity 
when measuring the shear profiles. Despite this simplification, 
our method is very satisfactory as we obtain a calibration bias 
(to) ^ -0.1 + 0.02 and a PSF residual cTc - 2.10"^, a value con- 



sistent with shot noise (see eq. 1 1 of Heymans et al. 2006 1. These 
two values can be compared to the results from other weak lens- 
ing methods that are summarized on Figure 3 of Heymans et al. 
(2006). So in the present case we introduce an underestimation 
of the true shear around 10% but with no significant systematics. 
In the rest of the paper we will therefore increase all measures of 
the shear by a 10% factor before we use them to derive cluster 
masses. 



3.3. Shear radial profiles 

To quantify the shape of a galaxy, we use the complex ellip- 
ticity defined by Bonnet & Mellier ( 1995 1. It relates the tensor 
of the galaxy brightness second moments to its shape. Applying 
the lensing transformation between the source and image planes. 
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Fig. 1. PSF treatment applied to stars: on the top-left panel we show the stars ellipticity before the Im2shape deconvolution. The 
corresponding distribution in terms of ellipticity components (el,e2) is shown in the top-right panel. Bottom-left panel shows the 
stars after the deconvolution by the PSF field (see text) and the bottom-right panel the corresponding distribution of (el,e2), which is 
more uniformly distributed than before the deconvolution. In both left panels is shown in red the scale corrsponding to a semi-major 
axis of 1". 



one find the simple relation ( Seitz & Schneider] 1997 1: 



1 -s-'e 



for 1^1 >1 

I e* -g* 



(1) 



where * denotes the complex conjugate, e**' being the complex 
ellipticity in the source plan and e in the image plan, g is the 
reduced shear 
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1 -K 



(2) 



which is a non-linear function of the two lensing functions: the 
(complex) shear y and the convergence k which is related to the 
projected mass density. 

We assume spherical symmetry of the mass distribution and 
we consider that unlensed galaxies are randomly oriented on the 
sky plane. This is observationally confirmed since the distribu- 
tion of the field galax y ellipticities can be fitted by a gaussian 
M0,cr2 ^ Q 252) ( Tyson & Seitzer|l988HBrainerd et al.|l996 i. 
In the weak lensing approximation (i.e. k «i 1) we get an un- 
biased estimator of the reduced shear by averaging the shape 
of background galaxies in concentric annuli around the cluster 
center Spherical symmetry also implies that the averaged ra- 
dial component of the shear vanishes, i.e. (e^) - 0, while the 
tangential component of the lensed galaxies traces the reduced 
shear (e^) = g. We also assume that (e*'') = with an intrinsic 



dispersion cr^ ^ 0.25. The measure of this averaged tangential 
component (e^) at different distances provides a reduced shear 
profile g{r) which is used to recover the total mass of the deflec- 
tor 

We fix the center of the lens at the position of the bright- 
est galaxy (BCG). The non-parametric mass reconstructions of 
the cluster sample that will be presented in Paper 2, show sep- 
arations between the cluster mass center and the location of the 
BCG smaller than 30" for all clusters except the most disturbed 
ones or those with a low S/N mass reconstruction. Real shifts 
between the BCG and the mass center are expected, especially 
in the most massive dark matter haloes or in non-relaxed clusters 
(see e.g. Skibba et al.,201 l.for a discu ssion of the central galaxy 
paradigm). However, Dietrich et al. (20111 also showed with 
simulated data that for mass reconstructions done with ground- 
based data, shifts are observed with a median value of 1' be- 
tween the true and the reconstructed mass center So here we 
can not reliably disentangle artefatcs in the 2D mass reconstruc- 
tions from real shifts between the BCG and the mass center. 
Moreover, given the redshifts and mass of the clusters, small 
miss-centerings disturb only the central data points in the shear 
profile. As discussed in Mandelbaum et al. (2010), excluding the 
central parts of the profiles in the fitting procedure reduces the 
underestimation of the mass because of miss-centering. Basic 
simulations (lensing a typical catalog of sources by a cluster rep- 
resentative of the sample in terms of mass and redshift) show that 
our fitting methodology gives consistent masses within their er- 
ror bars up top miss-centering of 1'. Therefore we consider that 
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Fig. 2. Distribution of the area (left) and orientation (right) of the stars selected in the field of RX J2228. 5+2036. The top panels 
represent the stars before the PSF deconvolution by Im2shape while the bottom panels show the distribution after deconvolution, 
with an average size close to pointJike objects and a uniform orientation of the stars. 
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Shear profiles are built in non-overlapping logarithmic an- 
nuli with rom = l-2r,>,, in order to have similar S/N ratios in 
each bin. Thanks to the wide field images taken with Megacam, 
we are able to measure shear profiles up to 600 to 800" from 
the BCG. We limit the fit of the shear profile in the central re- 
gion to an inner radius set at 25" (which roughly corresponds 
to the Einstein radius of the clusters) except for 4 clusters pre- 
senting a shear signal that is too disturbed in the central part. 
In these cases, the limit is fixed to the radius where the signal 
becomes significantly positive. This inner limit reduces the im- 
pact of miss-centering and avoids complications due to the non- 
linearity of the reduced shear with mass density. It also avoids 
the noise induced by galaxies which shape is badly estimated be- 
cause they are too close to bright cluster galaxies (e.g. |Cypriano| 
let aI.|2004D . 



Within these limits, 10 to 15 data points are available in the 
Fig. 3. Distribution of the errors on the galaxies ellipticity esti- sheai" profile. In order to improve the statistical weight of the 
mator given by Im2shape (here for RX J2228. 5+2036). The min- results, we have chosen to consti'uct 25 shifted profiles by mov- 
imum of the histogram is reached for an uncertainty 0-^ ~ 0.25. ing the position of the lower bound of the first bin by 1/25 of its 
This value gives a natural cut to clean up the catalogs of galaxies width. This allows to reduce sampling effects in the first points of 
used to measure the shear profiles. the profile where the averaged tangential ellipticity can signifi- 

cantly change according to the way the galaxies are binned. Each 
of these 25 profiles are fitted using a standard minimization 
(the distribution of field galaxies ellipticity being gaussian, 
in our case the BCG is a good estimator of the actual mass center e.g. Brainerd et al. (1996), the likelihood follows gaussian 
and that this assumption does not have any consequence on the statistics). Doing so, we limit the correlation that would arise 
total mass determination. with overlapping bins. 
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To improve the shear estimator, we weighted the average 
value of the tangential component of the galaxies ellipticity with 



the inverse of the ellipticities distribution variance Uj - l/cr^ 



IIS 



(3) 



(T II is the quadratic sum of the errors (p-ex^cTei) given by 
Im2shape for each object and the intrinsic dispersion of the 
galaxies ellipticities cr,„, = 0.25. Note that we also cut our 
catalogs of galaxies by limiting the uncertainties returned by 
Im2shape to (cTgi, crf,2) < 0.25. This value corresponds to the av- 
erage minimum in the uncertainties distributions (Figure [3]l, the 
objects above this threshold being most likely false detections or 
very faint/small galaxies. Note also that thanks to our weighting 
scheme their removal does not have any consequences on the 
final mass estimates (see Cypriano et al.|2004| for details). 

To assess the uncertainties on the best fit parameters, we pro- 
ceed as follow. Assuming that the measures ^/ have a gaussian 
uncertainty cr, (although the distribution of the errors appears to 
be highly non-gaussian on Figure [3] ellipticities for each galaxy 
generated by Im2shape have a gaussian distribution), we gener- 
ated 1000 Monte Carlo draws of each of the 25 shifted shear 
profiles, and applied the procedure. We then obtain 25000 
values of the best fit parameters. The mode of the resulting dis- 
tribution gives the most likely solution along with the 1 , 2 and 3 
cr asymetric errors. 

3.4. Normalization of the shear profiles 

To limit the contamination of the sources catalogs we introduce 
several magnitude and color cuts which are discussed in details 
in Paper 2. In summary, we limit the magnitude range to 21- 
m_comp+0.5 in r' and remove all objects located within a color- 
band defined by the bright ellipticals red sequence up to = 23. 
This reduces the galaxy density by about 1/3 even if some con- 
tamination by bluer galaxies remains. To overcome this problem, 
we followed the method proposed by Hoekstra (2007) which as- 
sumes that the number density of cluster members simply de- 
creases with radius as 1/r. We adjust the normalization by fitting 
the profile of galaxies excess relative to the background level and 
we correct the measured shear according to the distance to the 
cluster center (Figure|4]i. As the richness of a cluster scales with 
its mass, this contamination is higher for more massive objects. 
So we apply the renormalization procedure for each cluster in- 
stead of using an average correction as done in Hoekstra (2007). 

Photometric redshift information is also necessary to correct 
for two effects which change the normalization of the shear pro- 
files. The first one is an assessment of the fraction of foreground 
galaxies in the catalogs and the second one is the determination 
of the average redshift distribution of the background sources. 
To get a reliable faint galaxies redshift distribution, we used 
the photometric redshifts of galaxies in reference fields, neglect- 
ing the cosmic variance between different fields ([Cypriano et al.| 



2004{|Hoekstra|2007[pmousin et al.|2009[|Oguri et al.|2009| l. 
Observations in the CFHTLS Deep fields obtained with the same 
instrument form the ideal data set for our study (Ga vazzi~&| 
Soucail 2007| l. They are much deeper than the present obser- 
vations so we can consider that their catalogs are complete up 
to our magnitude limit. We used the photometric redshifts deter- 
mined with the HyperZ software (Bo lzonella et al.|2000) . They 
were calibrated and validated by comparison with spectroscopic 
redshifts obtained in the VVDS Dl field ( |Le Fevre et al.|2005| l 




800 



Fig. 4. Density profile of the background galaxies used for the fit 
of the shear profile in the field of MS 1621.5-1-2640. Despite the 
removal of the galaxies within the red sequence, a substantial 
fraction of cluster members remains. The density profile is fit 
(dashed red line) to give the boosting factor applied to the shear 
profile. 



and in the Growth Deep survey ( Weiner et al.|2005) . Once our 
photometric selection criteria were applied on the CFHTLS cat- 
alogs, we obtained a redshift distribution that is supposed to be 
representative of the catalogs we use for the weak lensing anal- 
ysis. 

We use this distribution of photometric redshifts to deter- 
mine the strength of the shear signal via the geometrical fac- 
tor /3(zl,Zs) - D IS I Dos {Dls is the angular-diameter distance 
between the lens and the source and Dqs the distance to the 
source). Both the convergence k and the shear y can be re-written 
as Poekstra et al.|2000| ): 



r 



(4) 



with Ps - yS/y6°° and y6°° , y" , /c°° are the corresponding functions 
for a source at infinity. In the weak lensing approximation, g ~ j, 
the redshift dependence of the signal is easy to take into account 
as the average reduced shear simply equals its value taken with 
the average geometrical factor. 

For galaxy clusters at low redshift, it is straightforward to re- 
place (fi) by the j6 value of an average redshift Zs, usually set to 
I ( |Okabe & Ume tsu 2008 ; Ra dovich et al. |2008| ). In the present 
case, the geometrical factor varies significantly with the sources 
redshifts so we computed the factor (J3(z)) for the whole dis- 
tribution of field galaxies. We account for the contamination 
by foreground galaxies given our selection criteria by setting 
PiZpiwt < Zciuster) - (valucs in a range of 20% to 30% de- 
pending on the cluster). Note that such an approach is only valid 
in the weak lensing regime because of the non-linearity of the 
reduced shear with p (see e.g. Seitz & Schneider] 1997 Hoekstra 
|et aLpOCjO *). However, as with remove the central parts from the 
fits, we expect these ifiiz)) to be a good approximation. The av- 
eraged geometrical factors are given in Table |2] along with the 
effective redshift defined as PiZeff) - (J3(z)). 
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Table 2. Summary of the weak lensing analysis. Columns are (1) cluster name; (2-3) NFW masses; (4) SIS velocity dispersion; (5) 
Einstein radius for a source at Zs - 2; (6-7) radius corresponding to the NFW masses; (8) average geometrical factor corrected from 
dilution (see text); (9) effective redshift derived from the average geometrical factor 
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4. Weak lensing masses 

4. 1 . Mass models 

We restrict our analysis to the two mostly used parametric mass 
models: the singular isothermal sphere (SIS hereafter) and the 
NFW model ( [Navarro et"aL]|1995| ). The SIS mass model is the 
most simple one to describe a relaxed massive sphere with a 
constant and isotropic velocity dispersion cr,,. The mass density 
profile writes as 



p(r) = 



and simple expressions are derived for the lensing functions 
y{r) = K{r) - RRl'^r where the Einstein radius scales as Re oc 
jScr^,. This mass model is a poor description of the actual mass 
distribution of clusters. So the fitting results of this mass model 
will not be used to derive the total clusters mass in the rest of 
this work. However it allows easy comparisons with other mass 
estimators such the dynamical analysis of the galaxies members 
that are also characterized by a velocity dispersion cr,,. It also 
provide an estimate of the Einstein radius which can be com- 
pared to results from strong lensing analysis. This latter quantity 
requires an estimate of the geometrical factor to be converted 
from cTi,. Values given in Table|2]are derived using a source red- 
shift = 2 which leads to radii in the range 20" - 40" (note 
that these values are not directly estimated from the signal in the 
central part as done in strong lensing analysis but derived from 
the fit of the weak lensing shear profile). 

The NFW mass profile ( Navarro eraL][T997l [2004| was in- 
troduced to model the mass distribution of a cold dark matter 
halo formed in a gravitational scenario of structures formation. 
It involves a shape parameter, the concentration c, and a normal- 
ization parameter often expressed as the scale radius rj: 



p(r) = 



Pc 



(r/r,)(l + r/r,)2 



Pc is the critical mass density of the Universe at the redshift of 
the cluster and 6c is a dimensionless factor which is related to the 
contrast density of a virialized dark matter halo A -p/pc = 200 
by 

_ 200 
^'-^ ln(l+c)-^ 




Mgoo with c=4 

Fig. 5. Comparison between the total mass M200 (expressed in 
units of lO'^h^gMg) obtained by fitting the shear profile with 
a NFW model. In the first case the concentration parameter is 
fixed to c = 4 while in the second case it is left as a free pa- 
rameter varying over the range 2 < c < 15, increasing the mass 
uncertainties by 25% typically. 



We also define R200 = c r^, often called the virial radius, and the 
mass enclosed within this radius 



M200 



800 TT 



PcR 



200 



A = 200 roughly corresponds to the virialized part of clusters. 
X-ray studies more often make use of A = 500, a density contrast 
more easily reached by the X-ray detection. Given the above re- 
lations, it is straightforward to analytically switch masses and 
radii from one A value to another. 

With two free parameters, the NFW model provides more 
freedom to adjust the mass profile of galaxy clusters. However, 
there is a well-known degeneracy between the total mass M200 
and the concentration parameter c when fitting the shear pro- 
file in the weak lensing regime. This comes from a lack of in- 
formation on the mass distribution near the cluster center The 
observed c - M degeneracy shows a much steeper slope than 
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the theoretical one ( |Okabe et al.|2010a] l and only a combination 
of strong and weak lensing can raise it and provide useful con- 
straints on the concentration parameter In the present case, with 
no strong lensing modeling of all the clusters in the sample, we 
decided to fix the concentration parameter and only fit the total 
mass M200- Indeed for most of the clusters, the fit of both pa- 
rameters leads to unrealistic values of c, either larger than 20 or, 
in most cases, sticked to the lower limit of the fit c = 2. Since 
a low concentration is associated to a lower mass when fitting 
a given shear profile, this leads to the apparent bias observed 
Figure |5] which is only an artefact due to poor constraints in 
the central parts of the shear profiles. For a more secure mass 
estimator we therefore fixed the concentration parameter c to a 
canonical value for massive clusters, C200 = 4, and fit the mass 
profile with only one free parameter. 

It is interesting to note that the ratio between the NFW mass 
M20() and the SIS mass integrated within Roqo is close to one for 
the whole sample: {M^f^^ / M'^ (R200)} = 0.98 + 0.06. This re- 
flects the fact that in this study of high redshifts clusters of galax- 
ies, weak lensing data only provide constraints on the total mass 
at large scales. It cannot efficiently constrain the radial mass 
profile contrary to clusters at lower redshifts where some dis- 
agreements are found between the 2 mass models ( |Okabe et al.| 
|20I0a» . 

In summary, we have in hands for each cluster the results of the 
fit by a SIS profile (cr,, and Re) as well as those by a NFW profile 
(M200 and /?20o)- The mass and the radius taken at A = 500 are 
also computed for comparison with the results of the X-ray anal- 
ysis. All these values are given in Table |2] Note that because we 
are interested here in the calibration of scalings relation and the 
scatter of individual clusters around the best fits, we do not have 
included the systematic uncertainty due to projection effects of 
the large scale structures. This effect should be less than 10% for 
rich clusters at intermediate redshifts ( |Hoekstra|2001| l. We also 
did not include the error on the geometrical factor (JS). Its value 
varies in a range 3% - 6% depending on the cluster when the dif- 
ferent CFHTLS Deep fields are used to estimate it. Moreover, in 
the present case, these 2 quantities are second order effects and 
the main source of uncertainties remains the intrinsic ellipticity 
of the lensed galaxies (e.g.|Hoekstra|: 



4.2. Comparison of the X-ray and weak lensing masses 

Raw XMM-Newton data were processed as described in 
IPointecouteau et al.| (|2005|l and [Pratt et al.| ([2007). As for the 
REXCESS sample (ICroston et al. 2008) 1, surface brightness 
profiles were PSF-coiTected, deprojected and converted to 3D 
density profiles using the non parametiic method described in 
Croston et al. (2006 1. 2D temperature profiles were extracted 
from XMM-Newton spectroscopy as described in Pratt et al. 
(|20T0l. The 3D profiles were recovered assuming analytical 
functions taking into account the PSF mixing effects and weight- 
ing the contribution of temperatures in various rings. A Monte 
Carlo procedure was used to compute the eiTors ( .Arnaud et al. 
2010 1. The hydrostatic mass equation was then applied to the 
density and temperature profiles and to their logarithmic gradi- 
ents to derive the observed mass profile for each cluster. Again 
a Monte Carlo method was applied to determine the eiTors on 
each point, taking into account over-constraints on the 3D pro- 
file imposed by the use of a parametric model. Observed mass 
profiles were finally fitted with a NFW profile, constraining the 
shape and the normalization parameters, i.e. C500 and Msoo. 

Figure |6] shows the comparison between the weak lensing 
and the hydrostatic equilibrium masses derived from the X-ray 




^^500 ^ 



Fig. 6. Masses obtained with the weak lensing analysis (y-axis) 
versus the hydrostatic equilibrium masses obtained with the X- 
ray analysis (x-axis). All masses are given in units of IO'^/i^^qMq 
and are computed at the density contrast A = 500. 



data for the 1 1 clusters. The comparison is made at a density 
contrast A - 500 to avoid any extrapolation of the X-ray data. 
The weak lensing masses M200 are simply converted into M500 
assuming a concentration parameter C200 = 4. This means that 
we are not comparing weak lensing and hydrostatic masses eval- 
uated in the same physical radius. We are instead focusing here 
on a raw comparison between fitting results. 
Over the whole sample, we almost reach a Icr agreement be- 
tween the 2 mass estimators since 7 clusters (i.e. ~ 64%) 
have compatible masses within their Icr interval of individual 
uncertainties, with a small scatter and no strong systematics. 
Fig. It] With the 1 1 clusters, we obtain an average mass ratio 
(M^'^/Mf "> = 1.92 + 1.58. Removing MS 1241.5-hl710 that is 
almost 3cr away from this average ratio, we find {Mwl/Mx) = 
1.49 ± 0.75, a value much less scattered. We recover here the 
'classical' result that is a moderate excess for the weak lensing 
masses. 

Similar studies were previously done on clusters at lower 
redshifts ( [Zhang et al.||2058t [Mahdavi et aL][2008l [Zhang et al. 
1 20 1 0[ ), but it is the first time we show this kind of comparison at 
such a high redshift range. Unfortunately, the weak lensing mass 
profiles are not sufficiently constrained to allow a more detailed 
study of the vaiiation with radius of the mass ratio M^l/Mx- A 
combination with strong lensing information would be needed to 
overcome this difficulty. It will be explored in a forthcoming pa- 
per for the EXCPRES clusters for which strong lensing features 
are detected (4 out of the 11 clusters). 

For the 4 outliers identified in Figure |6j their individual 
properties will be presented in Paper 2. However it remains 
difficult to fully understand on a case by case why these 
clusters show such discrepancies. It certainly results from 
the mixed effect of intrinsic physical departures from hydro- 
static equilibrium and the limited quality of the lensing data. 
Some differences are also statistically expected in unbiased 
samples of clusters and can be characterized with numerical 
simulations ( Meneghetti et al.|[2010[ l. However, because of the 
limited size of the sample, we did not separate dynamically 
perturbed clusters from the relaxed one and try to highlight vari- 
ations in the normalization and scatter of the Mwl - Mx relation. 
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-0.5 0.5 

(V-Mx)/V 

Fig. 7. Relative difference distribution between the weak lensing 
masses and the X-ray masses for the 1 1 clusters of our sample. 
Masses are given at the density contrast A = 500. 

5. Mass-observable correlations 

5.1. Methodology 

Because both masses and observables have uncertainties with 
possible correlations (e.g. the richness depends on /?20o and thus 
on the total mass), a simple minimization cannot be used to 
calibrate scaling relations. Moreover, we expect to have an in- 
trinsic dispersion around the main scaling relation which needs 
to be included and evaluated. So we preferred to perform a lin- 
ear regression using the orthogonal BCES estimators {Bivariate 



Correlated Errors and intrinsic Scatter, [Akritas & Bershady 
1996[ l. This approach has already been used by several groups 



for the calibration of mass-observable scaling relations ( Morandi] 
|et al.,,2007, Pratt et al. 2009 ), expected to be characterized by 
power laws and so fittied by linear relations y - ax + b m the 
log-log plan. For each relation, both variables are normalized 
by a pivot close to their mean. This insures that the logarithmic 
slope and normalization of the relation are nearly independent 
parameters. 

Once the best fit of the correlation is obtained, the total dis- 
persion over the points 



1 



Ypi - axi - Bf 



equals the quadratic sum of the statistical dispersion ctj. 
the intrinsic one cr,„,, which thus can be expressed as: 



1 



N r 



N-2 



iyi - axi -B) - 



N-2 



N 



(5) 



and 



(6) 



For each relation, we also estimated the Spearman correla- 
tion coefficient p: 

n{n^ - 1) 

where di is the ranking difference of the i-th sorted mass and ob- 
servable vector, p is a ranking correlation that shows the mono- 
tonic degree of the relation between the two variables (p = -1 
for a strictly decreasing relation, p = +1 in the opposite case). 
Compared to the classical Pearson coefficient, the Spearman co- 
efficient does not depend on the slope of the correlation. Such 




2.2 



Fig. 8. Mass-richness scaling relation calibrated with the NFW 
lensing masses M2m. The black line shows the best BCES or- 
thogonal fit. The region delimited by the two dashed lines corre- 
sponds to the intrinsic dispersion of the galaxy clusters around 
the best fit. The shaded area gives the statistical Icr uncertainty 
given by the best fit parameters. 



a coefficient has been used in similar studies to quantify the 
strength of a correlation, e.g. Lin et al. (2004 1. 



5.2. Optical scaling laws 

First, we compared the weak lensing mass M20i) to the optical 
richness A'200- We define A^200 as number of galaxies within 7?2oo 
that belong to the cluster red sequence (the precise definition 
will be presented in Paper 2). We keep only the brightes t mem- 
bers with L < 0.4-L*. L* is derived from the values of illbert et al.l 



(2005 1 and according to the cluster redshifts. Richnesses are cor- 
rected from the background contamination but are not corrected 
from any luminosity cut nor any color selection. The aim is to 
remain consistent with previous studies that focus on the bright 
elliptical cluster galaxies. 

The BCES calibration of this M2()o - A^2{)o relation, given in Table 
[3] presents a correlation with a logarithmic slope a - 1.04+0.38 
(Fig.[8]l. Despite a rather large uncertainty, the slope is very close 
to 1, in agreement with the simplest model of structures for- 
mation ( |Kravtsov et al.|p004| l. The uncertainties of the BCES 
fit are dominated by the intrinsic dispersion (dashed lines) due 
to the limited size of the cluster sample and increased by one 
outlier, RX J2228. 5+2036. These results are coherent with pre- 

20071 



Johnston et al. 2007| |Reyes et al.| 2008 


Rykoff et al. 


2008 


Mandelbaum et al.|2008 Rozo et al.|2009 


. We also compared 



the X-ray masses with the cluster richnesses. This leads to a cor- 
relation with a slope a - 1.14 + 0.18. In addition, we also looked 
at a possible correlation between the mass and the richness eval- 
uated in a given physical radius. In practice, we used the richness 
determined within 1 Mpc (roughly equal to 0.5/?20o)- The BCES 
estimator gives a slope a - 1.60 + 0.92, a result with a very large 
uncertainty that is useless with the present data. 

To compare these different approaches, we computed the rel- 
ative error R - \{Mwl - Mpfo_„)/MwL\ averaged over the sample 
of clusters after removing the main outlier RX J2228. 5+2036. 
The M-A^200 relation leads toR = 0.27 while with the M-Ni Mpc 
relation we find R - 0.41. Even if the conversion between rich- 
ness and mass within a fixed aperture is less physical than the 
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use of the virial radius R200, we show here that it is still accept- 
able. Clearly the uncertainties in the scaling relations are domi- 
nated by the intrinsic dispersion more than by the choice of the 
working radius (see also |Andreon & Hurn 2010). The accurate 
determination of /?200 is a second order effect and the use of a 
simple fixed aperture for the measure of the cluster richness can 
already lead to a correct mass proxy. This is particularly true as 
our sample covers a restricted interval in both mass and redshift, 
so the fixed aperture is not far from a fixed density contrast. 

Second, we tested the relationship between mass and op- 
tical luminosity. The M - L correlation gives a roughly con- 
stant mass-to-light ratio as the BCES estimator returns a slope 
a = 0.95 + 0.37. Note that to be consistent with other similar 
works, we do not correct the total luminosities from incomplete- 
ness. Our catalogs are limited to 0.4L* and a correcting factor of 
~ 1 .6 should be included to get the total luminosity, assuming a 
standard Schechter luminosity function. This is done in Paper 2 
which discusses in details the M/L ratio of the clusters. 
The correlation between X-ray masses and optical luminosities 
gives a = 1 .20 + 0.20, a value also compatible with a constant 
mass-to-light ratio. Note that, as for the M - N relation, the cal- 
ibration with the X-ray masses gives a smaller intrinsic disper- 
sion. 

Previous studies tend to find a slope in the M — L rela- 
tion sUghtly larger than 1 ( Lin et al.|2004 Popesso et al.|2005 



Bardeau et al. 2007, Reyes et al.||2008| l This tendency is pre- 
dicted by some of the semi-analytical models of galaxy forma- 
tion, at least for massive clusters ( |Marinoni & Hudson|2002| l in 
which the star formation efficiency can be inhibited by the long 
cooling times of hot gas on large-mass scales. Merging processes 
and galaxies interactions can also affect the luminosity of galax- 
ies and are more efficient in massive clusters. But because of the 
small mass range covered by our cluster sample, it is difficult to 
draw firm conclusions on any possible variation of the M/L ra- 
tio with mass. It would be necessary to add galaxy groups in the 
sample, although accurate mass determinations are quite difficult 
for low mass structures. 

Moreover, the complex method to determine the value of 
^200 represents an issue to estimate L200- So we repeated the 
same procedure as described above and compared the masses 
with luminosities integrated inside a 1 Mpc radius. We also 
calculated the average relative errors on the mass. We obtain 
R = 0.47 for masses derived from the M - L\ Mpc relation while 
we get R - 0.33 when we use the M - L200 relation. Again, 
the intrinsic dispersion is the main source of uncertainties, the 
precise determination of 7?2oo being a second order effect. 

5.3. X-ray scaling laws 

We focus here only on scaling laws fitted with ffie weak lens- 
ing masses taken at A = 200. The calibration of the Lx - M 
relation has been the subject of many studies, mainly at low 
redshifts. Mass estimators are usually derived from X-ray data 
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2009 


Vikhlinin et al. 


2009 1. Several attempts to use weak lens- 



ing masses were also proposed more recently ( Hoekstra]|2007 
Bardeau et al |2007[ jRykoffeFal . 2008; Leauthaud et al.|2010 
Okabe et al.||2010b||. Most these studies converge towards a 



mation based only on gravitation and assuming a fixed density 
contrast to derive masses ( |Voit|2005] l. 

Several processes can affect the properties of a cluster, both gen- 
erating an intrinsic dispersion around the expected Lx - M rela- 
tion and a breaking of self-similarity. Most of previous studies 
point towards a slope that is slightly larger than the expected one, 
i.e. a - 1.6-1.8 instead of 4/3. The intrinsic dispersion is found 
to be quite large, up to almost 50% (Reiprich & Bohr inger 2002]l. 
However, it appears that when cool cores are removed, this in- 
trinsic dispersion is smaller ( Alle n|1998 Markevitch| 1998f | Voit| 
et al. 2002; Maughan 2007; Pratt et al.|2009| l. Indeed, flieX-ray 
luminosity is strongly affected by the physical mechanisms rul- 
ing the baryonic content of clusters (pre-heating, radiative cool- 
ing, feedback of supernovae and active galactic nuclei ...), and so 
is the Lx - M with respect to the simplest gravitational model. 
These processes are modeled in up-to-date numerical simula- 
tions that combine the gravitational evolution of dark matter 
structures with the hydrodynamical behavior of the intra-cluster 
gas ( |Borgani et al.|2004[[Ka"y|2004[[Borgani|2008| l. These sim- 
ulations also seem to favor a slope of the Lx — M relation larger 
than the canonical one. 

The results of the BCES estimator for the Lx - A^^q relation 
are presented in Table [3] Although our Lx measurements do not 
purely equate Lx,200 (see Table[T]i, our best fit relation is in good 
agreement with previous works as we find a slope a - 2.06 + 
0.44. Despite a strong correlation, we also observe a large scatter 
around the best fit with cr,„, - 0.19. 

These encouraging results should nonetheless be taken with 
caution as our sub-sample of the EXCPRES clusters was de- 
signed to keep only the most luminous clusters in the sample. 
It therefore covers a limited range in X-ray luminosity. In or- 
der to explore a larger range, we combined this sample with 
clusters from two studies: the COSMOS galaxy groups analyzed 
in L eauthaud et al.| ( [2010|l and the maxBCG galaxy groups and 
clusters presented in |Rykoff et al. ( 2008| l (Fig.|9]). Both are using 
stacks of objects with similar richness to improve the statistics 
of the shear profiles for low mass objects. Their X-ray luminosi- 
ties are also stacked within /?200 apertures, whereas our values of 
the luminosities are computed within the detection radius (see 
Table [T]). However, we checked that the ratio between R^g, and 
/?2{)o for our sample is ~ 0.9 on average. We therefore consider 
that the present values are close enough to the values within 
R21X), thus allowing us to compare and combine them with the 
aforementioned works. This combination of the three data sets 
allows to cover an extended range in mass and luminosity. We 
fit them simultaneously and the best fit relation gives a slope of 
a - 1.57 + 0.07. This slope is lower than that from our data 
only but still in disagreement wiffi ffie hierarchical prediction of 
4/3. The intrinsic dispersion stays roughly unchanged over the 
combined sample. 

Although it is not nowadays much used as a mass proxy, 
we then checked the Mwl - T relation. We derived a slope of 
= 1.33 ± 0.42, a value compatible with the self-similar pre- 
diction of 1.5 within the Icr limit. This weak lensing calibration 
is compatible with those by Hoekstra ( 2007| l and Okabe et al 



power law ffiat differs from the relation expected in the hierar- 
chical scenario of structures formation where the X-ray luminos- 
ity scales with the mass as FZ^Lx oc {F~M,o,)*^^ . The F, factor 
stands for the evolution of the scaling relation and reduces to 
h{z), the reduced Hubble constant, for models of structures for- 



( 2010b| l, who respectively derive a slope of a = 1 .34 ± 0.29 and 
1.49 ±0.58. 

We also computed the average relative error R between the 
weak lensing masses and the mass proxies determined with X- 
ray data. The Lx -M relation gives R = 0.28 and the M -T lead 
to R = 0.31. These two values are of the same order of those 
previously obtained for the M - N and M - L correlations with 
optical data. 
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Table 3. Summary of the fitting results for the mass-observable scaling relations. Columns are (1) scaling relation; (2) best fit loga- 
rithmic slope; (3) best fit normalization; (4-5-6) total-statistic-intrinsic logarithmic dipsersions; (7) Spearman coefficient. Masses are 
expressed in units of 10'^ /i^^q Mq (except for the Lx- M relation with masses in units of 5 x lO'** /i^y Mq). The optical luminosities 
are given in units of 7 x 10'^ /i^g Lq, the richness is normalized to 80 and the temperature to 5 keV. X-ray luminosities are given in 
units of lO'*^ /ZyQ erg s"'. 



scaling law 




a 


A 






^stat 


a-int 


P 


- iVzoo 

^200 ~ -'^200 


1.04 
0.95 


±0.38 
±0.37 


1.27 ± 

1.28 ± 


0.12 
0.20 


0.17 ± 
0.18 ± 


0.05 
0.04 


0.11 ± 
0.10 ± 


0.01 
0.01 


0.13 ± 
0.14 ± 


0.06 
0.05 


0.66 
0.64 




2.06 
1.33 


±0.44 
±0.42 


0.29 ± 
1.03 ± 


0.07 
0.20 


0.27 ± 
0.15 ± 


0.06 
0.04 


0.16 ± 
0.09 ± 


0.03 
0.01 


0.19 ± 
0.12 ± 


0.08 
0.05 


0.85 
0.73 
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Rykoff et al. 08 O 
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Fig. 9. Weak lensing masses versus X-ray luminosities using the EXCPRES clusters (black boxes) combined with the stacked 
clusters and groups of Ryk off et al.| (2008) (green circles) and of [Leauthaud et aL] ( |2010| l (red crosses). The solid lines show the 
best BCES fit and the corresponding intrinsic dispersion for the three data sets. For comparison the BCES best fit obtained with the 
EXCPRES weak lensing sub-sample only is plotted (dashed line). 



6. Summary and conclusions 

In this paper we provide for the first time a thorough weak lens- 
ing analysis of a sample of 11 galaxy clusters (a sub-sample of 
the EXCPRES sample - Arnaud et al. in preparation) in a rela- 
tively high redshift range, 0.4 < z < 0.6. We have conducted a 
careful validation of our weak lensing pipeline, thanks in partic- 
ular to the STEP simulations which allowed to re-calibrate our 
shear measurements. The shear profiles were fitted with a NEW 
profile with a fixed concentration parameter C200 - 4. This is 



due to the lack of good constraints on the shear profile near the 
center so only the total mass M200 is a reliable measure. These 
weak lensing masses were compared to the X-ray masses de- 
rived from XMM-A^ewro« data. Over the 11 clusters, we find 
a good match between these two mass estimators for 7 clusters 
while we found 4 outliers, thus having almost a Icr agreement in 
the two methods. Such results are consistent with other similar 
studies performed for lower redshift samples. 
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The primary goal of this study was to look for correlations 
between the total mass of galaxy clusters and several baryonic 
tracers. In our sample of massive clusters, we found that the to- 
tal mass is proportional to the optical richness of galaxies. This 
result is in good agreement with previous works at lower red- 



shifts (|Becker et al. |2007| 


Johnston et al. 


2007 Reyes et al. 


2008VRykoff et al.||2008||M 


[andelbaum et a 


.||2008| |Rozo et al. 



2009 ). We reach the same conclusion for the correlation between 
mass and optical luminosity, with a constant ratio M/L, even if 
the correlation is weaker with a larger intrinsic dispersion. Using 
XMM-Newton data and the weak lensing masses, we build scal- 
ing relations with the X-ray total luminosity and the global spec- 
troscopic temperature. Both correlations, although weakly con- 
strained due to limited coverage of the mass range, give results 
in good agreement with previous studies: (i) the mass-luminosity 
relation presents a non self-similar slope, larger than expected 
from a purely gravitational model of structure formation, (ii) the 
mass-temperature relation is roughly compatible with the hier- 
archical prediction. We have extended our investigation of the 
Lx - M relation by combining our cluster sample with two sam- 
ples of groups and clusters for which weak lensing masses were 
obtained (Rykoff et al. 2008; Leautha udet al.|2010| l. The result- 
ing Lx - M2()o relation covers nearly two decades in mass and 
displays a regular shape with a well constrained slope. 

Even if the limited size of our sample and its reduced mass 
range coverage allow us to only probe the high end of the mass 
distribution of clusters, we have demonstrated the synergy with 
the X-ray measurements (i.e., mass, luminosity and tempera- 
ture). More specifically, we have confirmed that the scaling rela- 
tions with the lensing mass proxy are aheady in place at z ~ 0.5, 
with, at first order, no significant departures from the same rela- 
tions at lower redshift (assuming a self-similar evolution). The 
limited quality of the lensing data for our sample restricts the 
added value of a joint analysis, but the weak lensing and X-ray 
coverage of cluster prove to be complementary. Indeed X-ray 
data allow, to efficiently probe the central regions out to about 
^500- Beyond this radius the X-ray brightness rapidly drops, 
whereas the weak-lensing signal picks up to characterize the 
outer parts of massive halos. A perspective work would be to 
investigate the strong lensing signal (see Table[T]l of each cluster 
in our sample, in order to peiform a full lensing analysis from 
the inner part (with the strong lensing signal) to the outer parts 
(with the weak lensing signal). This would provide a coherent 
mass proxy estimator to be compared directly with X-ray and 
dynamical proxies. 

Finally, in the perspective of large optical surveys such as 
the Euclid space mission, self calibrated mass proxies will be 
needed for the full scientific exploitation of the mission. The M- 
N is an obvious candidate, which tight calibration will require 
validation against other relations and mass proxies, such as the 
X-ray mass. The M - N relation we provide in this paper is an 
early result in the redshift range 0.4 < z < 0.6. It illustrates the 
difficulty to constrain this particular scaling relation. Provided 
the estimator is made as universal as possible, which can be 
a non trivial task, this relation is promising in the framework of 
optical surveys. 
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Annexes 



(7^ = 1493 km.s 




Fig. 10. Reduced shear profile normalized by the virial radius 
g(r/R2Qo) for the galaxy cluster MS 0015.9+1609. The veloc- 
ity dispersion cr,, and the total mass Mioo fitted to this profile 
are given (red curve for the SIS model, blue one for the NFW 
model). 
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Fig. 11. Same as Figure 
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for the cluster MS 0451.6-0305. 
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Fig. 12. Same as Figure 



10 



0.5 1 2 

for the cluster RXC J0856. 1+3756. 
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Fig. 13. Same as Figure 



10 



0.5 

for the cluster RX J0943.0+4659. 




Fig. 14. Same as Figure 10 for the cluster RXC J 1003. 0+3254. 
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Fig. 15. Same as Figure 10 for the cluster RX Jl 120.1+4318. Fig- 18- Same as Figure[10|for the cluster RX J1347.5-1 144 
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Fig. 16. Same as Figure 
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for the cluster RXC J 1206.2-0848. Fig- 19. Same as Figure 
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for the cluster MS 1621.5+2640. 
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Fig. 17. Same as Figure 
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for the cluster MS 1241.5+1710. 
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Fig. 20. Same as Figure 10 for the cluster RX J2228.5+2036. 
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